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Using the Multislit Imaging Spectroscopy (MSIS) technique at the F0RS2 spectrograph on VLT-UTl, we have identified 
60 Intracluster Planetary Nebula (ICPN) candidates associated with the Intracluster Light (ICL) in the central region of 
the Hydra I cluster. Hydra I is a medium compact, relatively near (~ 50 Mpc), rich cluster in the southern hemisphere. 
Here we describe the criteria used to select emission sources and present the evidence for these PN candidates to be 
associated with the ICL in the Hydra I cluster. We also show, using the luminosity-specific planetary nebulae number, the 
a parameter, that the expected number of PNs associated with the stellar population of the central cD galaxy NGC 3311 
in the cluster is close to the number of PNs detected. 
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1 Introduction 

Diffuse intracluster light (ICL) has been observed both in 
nearby and intermediate redshift clusters (Mihos et al. 2005, 
Feldmeier et al. 2004). Observations show that the diffuse 
light in galaxy clusters amounts to ^ 10% of the total light 
emitted by the cluster galaxies (Aguerri et al. 2005; Zibetti 
et al. 2005). 

Cosmological simulations of structure formation predict 
intracluster stars to be lost from galaxies in interactions dur- 
ing the assembly of galaxy clusters (Napolitano et al. 2003; 
Murante et al. 2004; Willman et al. 2004; Sommer-Larsen, 
Romeo & Portinari 2005). Simulations suggest that most 
of the diffuse light in galaxy cluster cores originates from 
mergers that lead to the formation of the brightest cluster 
galaxy and of the other massive galaxies, while the tidal 
stripping mechanism dominates the formation of the ICL 
at larger radii (Murante et al. 2007). 

The study of the properties of this diffuse component 
has then an important role in understanding the mechanisms 
relevant for the evolution of galaxies in high density envi- 
ronments, and the formation history, dynamics and merging 
status of clusters. 

Wide field imaging has already shown the morphologi- 
cal complexity of the ICL: studies of the Virgo and Coma 
cluster ICL have shown that it is characterized by tidal fea- 
tures like streamers and extended galaxy halos (Adami et 
al. 2005; Mihos et al. 2005). Due to its intrinsically low 
surface brightness, /iy > 28.5 mag arcsec^^ (Feldmeier 
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et al. 2004; Mihos et al. 2005), the kinematics of ICL can 
only be studied using the Intracluster Planetary Nebulae 
(ICPNs) associated with this stellar component (Arnaboldi 
et al. 2004, Gerhard et al. 2007, Doherty et al. 2008). ICPNs 
are relatively easy to detect because their spectra are char- 
acterized by two strong emission lines: [OIII]A4959A and 
[OIII]A5007A, with relative flux ratio 1 : 3. 

The goal of our project is to measure the velocity dis- 
tribution of the PNs associated with the diffuse light in the 
central region of the nearby Hydra I cluster (Abell 1060, 
D ^ 50 Mpc, z ^ 0.0126), a medium compact, non-cooling 
flow, rich cluster in the southern hemisphere. This will be 
presented in a forthcoming paper. Here we present the se- 
lection criteria and the evidence for the detection of PNs 
associated with the ICL in Hydra I. 

At the distance of the Hydra I cluster, the flux of the 
0[III] A5GG7A emission line of a PN is less than 8 x lO^i* 
erg s^^cm^^, therefore we need to reduce substantially the 
noise from the night sky in order to detect these lines. This is 
possible by using the Multislit Imaging Spectroscopy tech- 
nique (MSIS: Gerhard et al. 2005, Arnaboldi et al. 2007). 
In what follows we define msoo? = — 2.5 log -Fsoo? — 13.74 
(the Jacoby formula, Jacoby et al. 1989), where F5007 is the 
integrated flux in the [QUI] A5007A emission line, and we 
assume a distance of 50 Mpc for the Hydra I cluster imply- 
ing 1" - 0.24 Kpc. 

In the next section we present the MSIS observations 
caiTied out with F0RS2 on UTl . In Section 3 we summarize 
the data reduction. We present the adopted selection criteria 
to identify the emission sources and the evidence for PNs 
associated with the ICL in the Hydra I cluster in Section 4. 
In Section 5 we show that the number density of PN candi- 
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dates detected is consistent with that expected for the stellar 
population associated with the central cD galaxy of the clus- 
ter, NGC 33 11. In Section 6 we summarize our results. 

2 Observations 

2.1 The Multislit Imaging Spectroscopy technique 

The MSIS technique consists of the combined use of a mask 
of parallel slits, a narrowband filter centered around the red- 
shifted [OIII] A5007A line, and a dispersing element. It is 
a blind search technique and allows one to obtain spectra 
of all PNs (and other emission line objects) that happen to 
lie behind the mask slits. Because the [OIII] emission lines 
from PNs are only ^ 30 km/s wide (Arnaboldi et al. 2008), 
when dispersed at spectral resolution R = A/ A A ^ 6000, 
their entire flux falls into a small number of pixels in the 
two-dimensional spectrum, determined by the slit width and 
seeing. By dispersing the sky noise on many pixels, the tech- 
nique allows measurement of very faint fluxes. We can de- 
tect PNs with magnitudes ~ 1.2 mag below the bright cut- 
off of the Planetary Nebulae Luminosity Function (PNLF; 
Ciardullo et al. 1998) and their positions and radial veloci- 
ties can be measured at the same time. 

2.2 Observational set up 

Data were collected in visitor mode during 2006 March 26- 
28, using the FORS2 spectrograph on UTl. The observed 
area covers the central region of the Hydra I cluster, around 
NGC 3311, at a = 10h36m42.8s, S= -27d31m42s (J2000). 

The FORS2 field of view (FOV) is 6'.8 x 6'. 8 wide, 
corresponding to ^ 10000 kpc-^, and is imaged onto a mo- 
saic of two CCDs, rebinned 2 x 2 in the readout. We used 
two narrow band filters, one centered at A = 5045A and 
a second one at A = 5095A, both with a FWHM of 60A. 
We cover, in this way, the whole range of cluster line-of- 
sight (LOS) velocities and, for fast, wlos > 4000 km s^^, 
and bright PNs, msooy < 29.3, the [OIII] A5007A emis- 
sion line falls in the red filter and the [OIII] A4959A line in 
the blue filter. Spectra were obtained with the GRIS-600B 
grism, which has a spectral resolution of 0.75 A/pixel at 
5075 A. The MSIS mask is made of 24 x 21 slits, each of 
them 0" .8 wide and 17" .5 long. Each slit is projected along 
the dispersion axis onto ^ 40 rebinned pixels. The effective 
area imaged by the slits is ~ 7056 arcsec^, that is ~ 4.5% 
of the whole FORS2 FOV. In order to cover the whole field 
the MSIS mask was stepped 15 times on the sky to cover 
the region between two adjacent slits. For each mask posi- 
tion 3 exposures, of 1200 sec each, were taken, ensuring a 
proper cosmic ray subtraction. The seeing during the three 
observing nights was in the range from 0" .6 to 1" .5. 

Arclamp calibration frames were acquired for the wave- 
length calibration, as well as flats and bias images. The flux 
calibration was done using long slit observations of the stan- 
dard star LTT7379 with narrow band filter plus Grism. 




Fig. 1 Upper panel: 2-D spectrum of a PN candidate. The 
dispersion axis is along the vertical direction. The emission 
line falls onto about 6 pixels both in the spatial and spec- 
tral direction. Lower panel: 1-D spectrum of the same PN 
candidate. The emission line has the same FWHM as the 
arc lamp lines, showing that the observed emission line is 
unresolved in wavelength. 

3 Data reduction 

Data reduction was carried out following the procedure de- 
scribed in Arnaboldi et al. (2007). After bias subtraction, 
the images were properly co-added and the continuum light 
from the two main galaxies was subtracted, with an fmedian 
filtering using a window of 19 x 35 pixels. 

Then we extracted and rectified the 2D spectra (Fig.[T]) 
of the emission sources and performed the wavelength and 
flux calibration. The total number of emission sources de- 
tected is 95. On the basis of the flux calibration the la limit 
on the continuum is ~ 7 x lO"^*' erg cm^^ s^^ A^^. 

A first classification of the emission line objects can be 
carried out according to the following criteria: 

- PN candidates: unresolved emission line objects, both 
in space and in wavelength, with no continuum; 

- background objects candidates: continuum sources with 
unresolved or resolved emission line. 

Of the 95 emission lines sources identified, 60 were clas- 
sified as possible PN candidates, the remaining sources as 
background object candidates. Monochromatic point like 
emissions appear in the final images as unresolved sources 
with a width of ^ 6 pixels both in the spatial and in the 
wavelength direction. The final spectral resolution is ^ 4.5A 
or 140 km s^^. Magnitudes for the PN candidates were 
computed using the Jacoby formula. 

4 Selection Criteria for the PN population in 
the Hydra cluster: the wavelength-magnitude 
plane 

We now describe the physical properties of PNs belong- 
ing to the Hydra I ICL in a two dimensional wavelength- 
magnitude space, as shown in Fig|2l 
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Fig. 2 The plot shows the wavelength vs magnitude space 
in which we describe the physical properties of a PN pop- 
ulation at 50 Mpc and bound to the Hydra I cluster Blue 
and red crosses are the PN candidates detected in the MSIS 
images. The various lines and the characteristics of the 4 
regions identified in the plot are explained in the text. 

A PN population is characterized by a bright cut-off 
of the PNLF, which according to Ciardullo et al. (1998) 
has absolute mag nitude of the [OIII] A5007A emission line 
of M* = -4.48. At the Hydra I cluster distance ( m ~ 
M = 33.49), this corresponds to an apparent magnitude 
of 7715007 — 29.1 (plotted as the dotted horizontal line in 
Fig|2|i. The apparent magnitude for the [OIII] A4959A emis- 
sion line is 1.2 mag fainter. 

The red inclined line in Fig |2] shows the apparent mag- 
nitude of a PN at the PNLF bright cutoff, as a function of 
the Hubble flow distance, and the blue line shows the same 
dependence for the [OIII] A4959A emission line. If the line- 
of-sight velocity corresponded only to the Hubble flow, then 
the magnitude of the two emission lines, [OIII] A4959A and 
[OIII] A5007A of a PN fainter than the bright cut-off would 
fall on the left side of the blue and red line, respectively. 

The systemic velocity of the Hydra I cluster f Hydra i = 
3683 km s^^ and its velocity dispersion (THydrai = 724 km 
s~^ (Christlein & Zabludoff 2003) identify a wavelength 
range in which [OIII] emissions from objects orbiting in 
the Hydra I cluster potential can be observed. At the sys- 
temic velocity of the cluster, the [OIII] emission lines are 
redshifted to A = 5068A and A = 5020A, respectively. 
In Fig|2]the central green line is at A = 5068A where the 
[OIII] A5007A line is redshifted to the cluster systemic ve- 
locity. The other two vertical green lines are at the wave- 
lengths bounding the ±3(t velocity range, respectively at 
A = 5032A and A = 5104A. 

The black continuous lines in Fig|2]show the measured 
profiles of the narrow band filters, T(A), normalized so that 
the maximum transmission is at the PNLF bright cutoff. 



They are centered respectively at A = 5045A and A = 
5095A and have a FWHM of - 60A. 

Considering the magnitude limit and wavelength range 
specified, we can identify different areas in this wavelength- 
magnitude space, separated by black horizontal lines: 

- Region #1 According to their velocity and magnitude 
all emission lines are [OIII] A5007A. 

- Region #2 It is the region where, considering the flux 
and wavelength, we can find both faint [OIII] A5007A 
or [OIII] A4959A emission lines. 

- Region #3 Here we have high flux and red wavelengths 
so that we can see the bright [OIII] A5007A emissions. 
For such emissions we expect to detect the correspond- 
ing [OIII] A4959A in region #2. 

- Region #4 Here the emissions are most likely [OIII] 
A5007A. In principle, we may find bofli [OIII]A4959A 
and [OIII]A5007A; however, if an emission were identi- 
fied as [OIII] A4959A its LOS velocity would be about 
8600 km s^^ which is more than 6ct away from the clus- 
ter systemic velocity. This PN would then not be bound 
to Hydra I and its velocity driven by the Hubble flow. 
Then its [OIII] A4959A magnitude should fall on the left 
side of the blue fine. Therefore [OIII] A4959A emission 
lines in this region are ruled out. 

We can now populate the wavelength-magnitude plane 
with the PN candidates detected in the MSIS spectra. This 
gives us important information to validate the PNs cata- 
logue. In the plot the blue and red crosses are the PN candi- 
dates detected in the wavelength ranges covered by the blue 
and red filters respectively. 

The first result is that the fluxes of these candidates are 
all consistent with those of PNs at the distance of the Hydra 
I cluster, in the range from 1.7 x 10^^® to 8.4 x 10"^^ erg 
s^^cm^^^ i 30 J > 7775007 > 28.9. 

The second is that for the two bright [OIII] A5007A 
emission sources identified in region #3, the corresponding 
[OIII] A4959A have been identified in region #2: these ob- 
jects are encircled in magenta. All others emission sources 
in the plot are [OIII] A5007A if in the cluster 

Moreover, we know that the la continuum upper limit 
flux in both filters is 7 x 10~ erg cm" -2 s-^A-i. Con- 
sidering that the compact HII regions detected in the Virgo 
cluster (Gerhard et al. 2002) have a V-band continuum flux 
of ^ 8.2 x 10~^^erg cm~^ s~^ A~^, corresponding to ^ 
9.6 X 10~^°erg cm~^s~^A~^ at the distance of the Hydra 
I cluster, this allows us to rule out from our sample com- 
pact HII regions such as or brighter than those observed in 
Virgo. Also, from the continuum upper limit flux we calcu- 
late that the Equivalent Width (EW) of the most luminous 
candidates is EW > 90A. This value is larger than those 
of [Oil] emitting background galaxies at z = 0.37, which 
have EW[oii] < 50 (Hogg et al. 1998). Therefore, we can 
rule out all contaminations to our sample except a few back- 
ground Lya galaxies. 

In summary, from analyzing the PN candidates in the 
wavelength-magnitude plane we learn that: 
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- fluxes are consistent with PNs at the distance of the Hy- 
dra I cluster; the only possible contaminants are Lya 
emitters at high redshift; 

- objects in magenta circles can reliably be classified as 
PNs because we are able to see both their [0111] A5007A 
and [OIII] A4959A emission lines; 

- for all other PNs we have detected the [0111] A5007A 
at the LOS velocity expected for objects bound to the 
Hydra I cluster. 

5 The PN population associated with the 
stellar population of the Hydra I cluster 

We now evaluate the expected number of PNs associated 
with the total stellar luminosity of NGC 3311 and we test 
whether the number of detected candidates is consistent with 
that expected given the sampled light, the flux limit and 
the average color of the stellar population. The luminosity- 
specific planetary nebulae number, the so called a parame- 
ter, specifies the number of PNs associated with the amount 
of light emitted by the stellar population of a galaxy (Jacoby 
et al. 1980). According to Buzzoni, Arnaboldi, & Corradi 
(2006) this is defined as a = Npm/Lb, where Np^ is the 
number of PNs of a galaxy and is is its total stellar lumi- 
nosity in the B-band. As shown in Buzzoni et al. (2006), the 
a parameter is correlated with the {B — F) colour index of 
the galaxies. 

To compare the theoretical number of PNs expected for 
NGC 33 1 1 in Hydra I wifli flie number detected in flie MSIS 
images we have gone through the following steps. We took 
{B-V) = 1.04 for NGC 33 1 1 from Prugniel & Heraudeau 
(1998) and flie colour excess E(B -V) = 0.079 from the 
NED database to determine the [B — V^q = 0.961 colour 
corrected for extinction. From this value we obtained a = 
1.5 X 10"'^ using Table #6 in Buzzoni et al. (2006). To calcu- 
late the predicted number of PNs in the MSIS images a com- 
pleteness correction is necessary. Since with our technique 
we detect only PNs down to 1.2 magnitude from the PNLF 
bright cut-off, the number we need is ai.2- This is simply a 
fraction of the total luminosity-specific PN number: ai 2 = 
a/30 = 5 X 10-9 (Buzzoni et al. 2006). The total B-band 
luminosity of the diffuse light from NGC 33 11 's envelope in 
the FORS2 field, excluding the central 34" where the galaxy 
is brighter than the sky, is Lb = 6.8 x 10^°^© (Mb = 
-21.7 from Vaesterberg, Lindblad & Jorsater 1991). Then 
the total number of PNs expected in the FOV is: Np^ = 
Oii.2/LB ~ 114. Because each MSIS mask covers an area 
of about 4.5% of the total FORS2 area, we thus expect to 
see 114 x 0.045 ^ 5 PNs per mask, i.e. ~ 75 PNs in the to- 
tal MSIS image sample. This is consistent within 1.5ct with 
the actual number of PN candidates, 60. 

6 Conclusions 

In this work, we have presented MSIS observations of PN 
candidates associated with the diffuse light around the cen- 



tral cD galaxy NGC 3311 in the Hydra 1 cluster. We have 
discussed criteria used for selecting these emission sources 
and have analyzed their properties in the wavelength-magni- 
tude plane. In total, we have identified 60 PN candidates 
around NGC 331 1, consistent with the number of PNs ex- 
pected from the diffuse envelope of this galaxy. 

In the next steps of our analysis we will focus on the 
kinematics of the diffuse stellar population in the halo of 
NGC 33 1 1 . We wiU study the histogram of tiie LOS veloc- 
ity distribution of the PNs candidates to obtain information 
about possible substructures in the ICL and the merging sta- 
tus of the cluster. Together with the X-ray temperature pro- 
file (Yamasaki et al. 2002) and the surface brightness pro- 
files of the stellar light this will give us constraints on the 
orbit distribution in the halo of this cD galaxy. 
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